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Abstract 

The phase space dynamics of dissipative quantum systems in strongly con- 
densed phase is considered. Based on the exact path integral approach it 
is shown that the Wigner transform of the reduced density matrix obeys a 
time evolution equation of Fokker-Planck type valid from high down to very 
low temperatures. The effect of quantum fluctuations is discussed and the 
accuracy of these findings is tested against exact data for a harmonic system. 
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Introduction-Dynamics in strongly condensed phase can be found almost everywhere 
in nature, e.g. for tunnel diodes in mesoscopic physics, for macromolecules in biological 
and soft matter systems, and for chemical reactions, to name but a few. Associated with 
this bunch of realizations is a rich phenomenology comprising prominent effects such as 
stochastic resonance |IJ, resonant activation 0, transport in ratchets and adiabatic 
electron transfer [Q which have been explored extensively in the last decade. However, most 
of these studies have focused on the domain of classical physics and less is known about 
low temperature quantum properties. This is mainly due to the fact that in contrast to the 
Fokker-Planck equation for the classical phase space distribution || , a simple time evolution 
equation for the reduced density matrix of a dissipative quantum system does not exist in 
general || . Typically, quantum fluctuations appear on a time scale h/3 ((3 = 1/ k B T) so that 
at lower temperatures the quantum stochastic process becomes strongly non-Markovian and 
intimately depends on the initial correlations between system and heat bath. An exact 
expression for the density matrix is available within the path integral formulation, but 
numerical evaluations are in most cases prohibitive ||. 

In the low friction limit, e.g. for quantum optical systems, progress can be made by a 
Born-Markov approximation leading to well-known master equations. The opposite limit of 
strong damping has only recently attracted considerable attention fTj- |10| . In particular, it 
was shown || that for strong friction and low temperatures the density matrix is essentially 
restricted to its diagonal part, the position probability distribution, and that this part obeys 
a quantum analogue of the classical Smoluchowski equation. The question we address here 
goes far beyond: Is there a time evolution equation in whole phase space valid from high 
down to low temperatures provided friction is sufficiently strong? The answer is yes and we 
will show this in detail below. 

Path integral formulation-The inclusion of dissipation within quantum mechanics is well 
established and is based on a system-l- reservoir formulation 0. The dynamics of the corre- 
sponding density matrix starting at t — from a general initial state W(0) reads 
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W{t) = exp(-iHt/h) W{0) exp(iHt/h) (1) 

where the Hamiltonian H = H$ + Hr + Hi contains a system, a reservoir (heat bath), 
and a system-reservoir interaction part, respectively. Dissipation arises as the effective 
impact of the reservoir degrees of freedom on the system dynamics within a reduced picture 
p(t) = trn{W(t)}. In the standard description the Gaussian statistics of the heat bath is 
modeled by a quasi-continuum of harmonic oscillators bilinearly coupled with the system. 
In fact, one this way regains in the classical limit the generalized Langevin equation of 
Brownian motion. The only non-perturbative treatment of the system-reservoir coupling is 
an exact elimination of the bath degrees of freedom by means of the path integral approach. 
For the position representation of the reduced density matrix one obtains 

p{Qf,q'f,t) = J dq.dq'i J{qf,q'f,t,qi,q^ A(fc,q£). (2) 

Here, the propagating function J(-) is a threefold path integral over the system degrees of 
freedom only. The two real time paths q(s) and q'(s) connect in time t the initial points 
qi and q\ with the fixed end points qf and q'f, while the imaginary time path q runs from 
qi to q[ in the interval %(3. The contribution of each path is weighted by exp(zS[g, q', q]/h) 
with an effective action S[g, q', q] not specified here explicitly. Basically, it comprises the 
actions of the bare system in imaginary and real time, respectively, and additional interaction 
contributions (influence functional), non-local in time, which in the reduced picture rule the 
influence of the reservoir onto the system. In the limit t — > one has </(<?/, q'f, t, q i: g-) — > 
PpiQi, Qi) S(q f - qi) 5(q' f - <?■) so that 

PiVf, <l'f, 0) = ppifif, q' f ) X(q f , q' f ) (3) 

with the reduced equilibrium density matrix P/3(q,q f ) = (q\tr R exp(—j3H)\q') and a prepara- 
tion function A characterizing initial deviations from thermal equilibrium. The crucial point 



here is that Eq. (H) is a correlated initial state [11|, in contrast to a factorizing initial state 



used in ordinary Feynman- Vernon theory. In the strong friction regime the latter one is 



not applicable in contrast to what has been done in |TC| . Even in the classical limit does a 



product initial state not allow for a consistent derivation of the well-known Smoluchowski 
equation. 

Overdamped quantum systems-The stochastic dynamics of dissipative quantum systems 
determined by Eqs. (0) and ([$]) is rather subtle as it strongly depends on initial correlations 
between system and reservoir as well as on the non-local self-interactions contained in the 



influence functional fly]. Both are governed by the damping kernel 



Km =r^i M ^f-j^ ( 4) 

Jo 7r smh{unp/2) 

where 9 = s — ir, < s < t, < t < h/3 and is the spectral density of the heat bath. 
Accordingly, the impact of the bath is completely specified by temperature T and I(uj). In 
particular, the spectral density is related to the macroscopic damping via 

= m / ~z — r cos ( ws ) • ( 5 ) 

Now, we turn to the simplifications that arise in the overdamped regime. For this purpose 
we define a typical damping strength as 

with j(ou) the Laplace transform of 7(i). For instance, in the ohmic case I(co) = M^fuj 
and also for the more realistic Drude model I(uS) = M^ujI/ {oj 2 + u 2 ) (cut-off frequency u c ) 
one finds 7 = 7. Given a typical frequency u of the bare system, e.g. the ground state 
frequency, by strong damping we then mean 

1 11 

J_»^,__. (7) 

^0 u c 7 

In other words, we assume the time scale separation well-known from the classical over- 
damped regime and extend it to the quantum range by incorporating the time scale for 
quantum fluctuations h/3. Correspondingly, we consider the dynamics Eq. (0) on the coarse 
grained time scale s 3> h/3, - and a 3> -. The consequences are substantial: (i) the 
strong friction suppresses non-diagonal elements of the reduced density matrix during the 



time evolution. This simply reflects the fact that a quantum system behaves more classically, 
the stronger coherences are destroyed by the presence of a heat bath, (ii) The real-time part 
K(s) of the damping kernel becomes local on the coarse grained time scale so that a time 
evolution equation of the form p(t) = £ p(t) with a time independent operator C may exist. 

Following the above simplifications the path integral formulation now allows for a pertur- 
bative treatment in the strong damping limit. To this end it is convenient to introduce sum 
and difference coordinates for the imaginary time path, f = (q+q')/2 and x = q—q', and sum 
and difference paths in real-time, r(s) = [q(s) +q'(s)]/2 and x(s) = q(s) —q'(s), respectively. 
The idea is to evaluate the path integrals in the sense of a semiclassical approximation by 
assuming self-consistently that non-diagonal elements, i.e. x and x(s) dependent terms, re- 
main small during the time evolution. Hence the effective action E[r, x, q] is expanded up to 
second order in the x coordinate of the imaginary time path and in the excursions x(s) of the 
real-time path integrals. Doing so we take sufficiently smooth potentials for granted. While 
this procedure is applicable to a wide range of spectral bath densities, we concentrate in the 
sequel on the quasi-ohmic case with a very large cut-off frequency u c ^> 7. It is worthwhile 
to note that we do not restrict the value of 7^,/? meaning that our analysis covers a broad 
temperature range from the classical (jfiP <C 1) to the deep quantum domain {p{Tif3 3> 1). 

Equilibrium density matrix-The question to what extent the thermodynamic equilibrium 
is affected by the strong damping limit, particularly at lower temperatures, is a very crucial 
one and also serves as a basis for the more involved treatment of the dynamical case. In 
terms of a sum over paths from q to q' in the time interval [0, H/3] the unnormalized reduced 
equilibrium density matrix in position representation reads 



with the bare Euclidian action Se[q] — Jo dr [Mq /2 + V(q)} of a particle of mass M 
in the potential V(q) and the Euclidian influence functional (f>E[q] = Jo dr J daKiir — 
ia) q{r) q(u). To solve the path integral we put q = g ma + Sq where q ma ,(cr) is the minimal 
action path to (H) and obeys the boundary conditions q(0) = q, q(h(3) = q', while 5q(o~) de- 
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(8) 



note quantum fluctuations with Sq(0) = 8q(%f3) = 0. The most convenient way to calculate 
g ma is to switch to Fourier space with respect to the Matsubara frequencies u n = n2ir/hf3. 
In the corresponding equation of motion to leading order the friction term prevails and one 
finds g ma = f (recall that x is assumed to be small, namely, of order I/-/7). Dynamical 
contributions, however, must be taken into account if we are interested in next order correc- 
tions. The final result for g ma including corrections of order I/7 is lengthy and thus omitted 
here; the minimal action takes the form 

S E (f, x)/h = pV(f) - A/3 2 V'(f) 2 + i 2 + 0{A/j) (9) 

with 

A = J_f I M 2M- 7 K) 

Apparently, A measures the typical strength of quantum fluctuations in position space, while 
Q is via the identification x/h — > p associated with the variance in momentum Q = (p 2 ). In 
case of Drude damping with a high frequency cut-off uo c both A and Q can be expressed in 
terms of \l/ functions. Then, for high temperatures ^U(3 C 1 we find A w h 2 f3/12M and w 
Mj f3. The friction dependence appears as a genuine quantum effect for lower temperatures 
and for jh(3 > 1 one has A w (h/Mjir) \og(-fhf3/2n) and Q w (MUj/'k) log(a; c /7). With 
increasing 7 the strong squeezing of quantum fluctuations in position induces enhanced 
fluctuations in the momentum (see fig. |l|), thus suppressing nondiagonal elements in the 
density matrix. 

By expanding the full action in (|8[) up to second order in Sq, the contribution of quantum 
fluctuations to (§) is obtained. Due to the strong friction higher order terms are negligible. 
The corresponding Gaussian integral is again calculated by utilizing a Fourier expansion in 
the Matsubara frequencies. Eventually, the equilibrium density matrix in the strong friction 
limit is found as 

Pf3 (x,f) = l e -mf)-ns 2 /27i 2 e Af3[f3V(f) 2 /2-3V"(f)/2] ^ 

z 
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where Z denotes a proper normalization factor, e.g. Z = J dqpp(0,q). Interestingly, the 
probability distribution is Gaussian in x, i.e. its Wigner transform (x/K — > p) Gaussian in 
momentum, even at low temperatures. Anharmonic corrections in x to the exponent are 
at most of order I/7 2 . The expression ( O ) is an important result since numerically exact 
calculations of the path integral (^) e.g. via Monte Carlo techniques become more expensive 
at lower temperatures. It further reveals that for strong friction the equilibrium density 
consists of a part which in phase space takes the form of a classical distribution, however, 
with an h dependent (p 2 ) and a part with A dependent quantum corrections (see also fig. [[]). 

Quantum Fokker-Planck equation-Fov the time evolution of the density we basically apply 
the same kind of semiclassical analysis as briefly outlined above. Since on the coarse grained 
time scale the dynamics is effectively Markovian it suffices to calculate the propagating 
function J(-) for a time step from t to t + St where St obeys 7/^0 ^> St h/3, 1/7, 1/uj c . 
Our goal is to derive from this result a time evolution equation for the density which after 
a Wigner transform gives rise to a quantum Fokker-Planck equation. 

We start by expanding the effective action up to second order in x and x(s). Specifically, 
for the potential terms in the real time actions one writes V(r + x/2) — V(r — x/2) = 
V'(r)x + 0(x 3 ). The complication that arises in evaluating the minimal action paths is that 
due to the correlations between system and bath the corresponding equations of motions are 
coupled. For anharmonic potential fields this necessitates in general a numerical evaluation. 
However, here it turns out that one may write g ma = r, + Sq and r ma = Tf + Sr where 
Sq and Sr are of order A or smaller. This way, Sq is determined similar as for the static 
case. By approximating V'(r) = V'(r/) + V"(rf) Sr + 0(Sr 2 ) the real-time paths Sr and x 
run effectively in a harmonic oscillator potential with frequency JV"(rf)/M subject to an 
external force V'(r/). The corresponding part of the propagating function has been derived 
in JO]] to which we refer for further details. Hence, to gain a time evolution equation for 
p(rf,Xf,t) in the form p{t) = Cp(t) we look for an operator C = C(xf,rf,d/dxf,d/drf,t) 
being at most second order in the coordinates and derivatives. Strong friction forces forward 
and backward real time paths to run very close to eachother so that quantum fluctuations — 
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responsible for diffusion — reduce to Gaussian noise. Below we will briefly discuss higher 
order corrections. To proceed, we make a general ansatz for £ with coefficients specified 
by comparing p with Cp. Eventually, we switch to classical phase space {xf,rj} — ► {p, q}, 
i.e. p(xf,rf,t) — > W(p, q, t), by the replacement Xf — > ihd/dp, r/ — > g and d/dxf — > 
d/drf — > d/dq. This leads to the main result of this article, namely, the time evolution 
equation for the phase space distribution £) of a dissipative quantum system in the 

strong friction limit (0) 

|^«.*)={|W.(«)+7ri-^ + 7< P »>J 

+ d^p [ 1 //3 + A^"(?)-(P 2 )/M]| W(p,g,f). (12) 

Here, we have introduced an effective potential V e s = V+(A/2) V", and A and also (p 2 ) = 
are specified in (fL0|). The first line on the r.h.s coincides with a classical Fokker-Planck 
operator in an effective force field [§], the second line describes quantum mechanical coupled 
p-q diffusion. In the high temperature limit ^hf3 — > the quantum Fokker-Planck equation 
(QFP) tends to the classical Kramers equation 0. For small but finite 7^/5 and in case 
of a harmonic potential the QFP coincides with the master equation gained by Haake and 
Reibold [T^j , but differs from the Caldeira-Leggett master equation |TB| by the p-q diffusion 



term [|Tj]. However, while these known master equations are restricted (for 7/cJo > 1) to the 
range of small 7^/?, the new QFP is valid for all ~fh/3. We only mention here, that the above 
QFP is not of Lindblad form due to the coarse graining procedure on which its derivation is 
based. Of course, the equilibrium solution to ( |12"D is given by the Wigner transform of (fLip . 
Let us briefly touch the question about higher order diffusion terms to (|T2|) . For harmonic 



systems they do not occur so that the QFP is in this sense exact |T5[ . In case of anharmonic 
potentials they result from non-Gaussian quantum fluctuations attributed to higher than 
second order derivatives in V(q). A rough estimate shows that anharmonic terms in Xf 
(leading to higher than second order derivatives in p) in the crucial low temperature range 
"fhp 3> 1 are of order l/[^7 3 / 2 log(u; c /7)] compared to the leading terms. 

Applications-With the QFP at hand, we are now able to study phase space properties 
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of quantum systems in strongly condensed phase. In particular, classically, due to the 
time scale separation the strong friction limit allows for a reduction of the phase space 
Kramers equation to the position space Smoluchowski equation p|,|T6f which has been of 
great importance in a variety of systems in physics and chemistry (cf. the Introduction). 
Its generalization to the low temperature quantum domain has been found recently in || 
by focusing on the path integral expression for the diagonal part p(q,q,t). Now, from 
the quantum phase space dynamics this latter result can be rederived and especially the 
influence of inertia effects can be explored. For this purpose we employ the projection 
operator techniques invoked in |16| to systematically reduce the QFP to position space. 
Along these lines we introduce the operators V = fp{p) J dp and Q = 1 — V where fpip) is 
the normalized momentum distribution in equilibrium according to fllTI). The next steps are 
straightforward and not presented here in detail. After some algebra one arrives to order 
A/7 3 at an equation for the position distribution n(q, t) = f^ 1 VW of the form 

" = W!T q { 1 + Wf W + Ar 'l } L ^ 4 » ' < 13 > 

Lqse = Kg + d/dq[l//3 + AV"} is the quantum Smoluchowksi operator already derived in 
H (slightly generalized to all values of 7^/3). A classical (A = 0) inertia correction oc V" 
appears, while quantum fluctuations enter through third and forth order derivatives of the 
potential. Overdamped quantum Brownian motion in position space thus becomes much 
more sensitive to the details of the potential profile at lower temperatures. 

There is one non-trivial case where analytical results are available and can be compared 
with exact ones, namely, the damped harmonic oscillator. As already mentioned above, in 
the quantum case initial correlations between system and reservoir are of crucial importance 



and render the calculation even for the harmonic case quite cumbersome [DJ. This is 
particularly true for the strong friction range where factorizing initial states cannot be used. 
We illustrate this in detail by considering the relaxation of expectation values in position 
(q(t)) and (q(t) 2 ) (see fig. |2|). Already for a moderate damping is the QFP solution for 
the mean position in very good agreement with the exact dynamics, while that based on a 
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factorizing initial state is completely off. Upon closer inspection one finds that in the latter 
case the g(0)-dependence is suppressed by a factor I/7 2 compared to the former ones due 
to the lack of initial correlations. For the mean square in position deviations between the 
QFP result and the exact one are almost invisible. Remarkably, the impact of quantum 
fluctuations due to A and (p 2 ), in the QFP encoded in the q-p diffusion term, is quite 
substantial. 

Conclusion-We have analyzed the time evolution of strongly damped quantum systems in 
phase space starting from the exact path integral expression in position space. A suppression 
of quantum fluctuations in position is associated with enhanced fluctuations in momentum. 
Our central result is a quantum Fokker-Planck equation applicable from the classical high 
temperature {"ffifl <C 1) to the low temperature quantum {^Tif3 ^> 1) range. This way, while 
in the weak friction limit master equations have been known for years, here, we have given 
the missing complement for strong damping. The QFP now opens the door to explore phase 
space features of many systems well studied in classical physics also at lower temperatures, 
e.g. Kramers rate theory, driven transport, or soft matter problems. 
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FIGURES 

FIG. 1. A as a measure for position fluctuations and = (p 2 ) vs. inverse temperature for 
different 7 and a Drude model with uj c /ujq = 50. A is scaled by H/Mujq, by KMujq. 

FIG. 2. Relaxation of position expectation values for a damped harmonic oscillator within 
a Drude model uj c /^o = 50, 7/^0 = 5, (jJq%(3 = 1. Initial values are for left panel: 
(p(0))/Mu {q(0)) = 1; for right panel: (<? 2 (0)) = (p(0) 2 )/M 2 u; 2 = (<? 2 )/4 (other moments are 
zero) where (q 2 ) is the quantum equilibrium variance. 
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